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Abstract. We continue our study of a magnetic Schrodinger operator on a 
two-dimensional compact Riemannian manifold in the case when the minimal 
value of the module of the magnetic field is strictly positive. We analyze the 
case when the magnetic field has degenerate magnetic wells. The main result 
of the paper is an asymptotics of the groundstate energy of the operator in the 
semiclassical limit. The upper bounds are improved in the case when we have 
a localization by a miniwell effect of lowest order. These results are applied 
to prove the existence of an arbitrary large number of spectral gaps in the 
semiclassical limit in the corresponding periodic setting. 



1. Preliminaries and main result 

Let M be a compact connected oriented manifold of dimension n > 2 (possibly 
with boundary). Let g be a Riemannian metric and B a real- valued closed 2-form 
on M. Assume that B is exact and choose a real- valued 1-form A on M such that 
e?A = B. Thus, one has a natural mapping 

mh> ihdu + Am 

from C£°(M) to the space Vi).(M) of smooth, compactly supported one-forms on 
M. The Riemannian metric allows to define scalar products in these spaces and 
consider the adjoint operator 

(ihd + A)* : f^(M) -> C C °°(M). 

A Schrodinger operator with magnetic potential A is defined by the formula 

H h = (ih d + A)*(ih d + A). 

Here h > is a semiclassical parameter, which is assumed to be small. If M has 
non-empty boundary, we will assume that the metric g and the magnetic potential 
A are smooth up to boundary and the operator H h satisfies the Dirichlct boundary 
conditions. 

We are interested in semiclassical asymptotics of the low-lying eigenvalues of the 
operator H h . This problem was studied in [4j [U UHl HU H21 [13 HH HH HU [20l EU 
[21 [11 [21 [25] (see 0[S] for surveys). 
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In this paper, we study the problem in a particular situation. We suppose that 
M is two-dimensional. Then we can write B = bdx g , where b £ C°°(M) and dx g 
is the Riemannian volume form. Let 

bo = min \b(x)\ . 

We furthermore assume that: 

• b > 0; 

• the set {x £ M : \b(x)\ = bo} is a smooth curve 7, which is contained in 
the interior of M; 

• there is a constant C > such that for all x in some neighborhood of 7 the 
estimates hold: 

(1.1) C- l d(x, 7 ) 2 < \b(x)\ -b < Cd(x, 7 ) 2 . 

The main purpose is to give an asymptotics of the groundstate energy of the 
operator H h . Denote by N the external unit normal vector to 7. Let N denote the 
natural extension of N to a smooth normalized vector field on M , whose integral 
curves starting from a point 1 in a tubular neighborhood of 7 are the minimal 
geodesies to 7. Consider the function fi 2 on 7 given by 

(1.2) p 2 (x)=N 2 \b(x)\, x£ 1 . 
By (|l.ip . it is easy to see that 

P 2 (x)>0, X £ "f . 

Theorem 1.1. 

There exists ho > 0, such that, for any h g]0, ho], 

X (H h ) = hb + h 2 ^+O(h 17 / s ). 
4&o 

where 

fio := inf j3 2 (x). 

The paper is organized as follows. In Section [2] we construct approximate eigen- 
functions of the operator H . This allows us to prove an upper bound in Theo- 
rem [O] (see Corollary I2.2|) . In Section [3] we prove a lower bound for Xo(H h ) and 
complete the proof of Theorem ll.il In Section|4l assuming the existence of a unique 
non-degenerate miniwell for the reduced spectral problem on 7, we construct more 
refined approximate eigenfunctions of the operator H , improving the upper bound 
of Section [2 In Section [S] we consider the case when the magnetic field is periodic. 
We combine the constructions of approximate eigenfunctions given in Sections [5] 
and [4] with the results of [7] to prove the existence of arbitrary large number of 
gaps in the spectrum of the periodic operator H h in the semiclassical limit. 

Acknowledgements. 
We wish to thank the Erwin Schrodinger Institute in Vienna and the organizers 
of the Conference "QMathll — Mathematical Results in Quantum Physics" in 
Hradec Kralove for their hospitality and support. 

2. Upper bounds 
2.1. Approximate eigenfunctions: the main result. 

The purpose of this section is to construct approximate eigenfunctions for the op- 
erator H h . Denote by R the scalar curvature of the Riemannian manifold (M,g). 
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Theorem 2.1. 

For any a; £ 7 and for any integer k > 0, there exist C and ho > 0, such that, for 
any h G]0, ho], there exists € C£°(M), 7^ , swc/i i/mf 

\\H h ^-\ h (k,x)t$\\<Ch 17 '*\\1%\\, 

where 

(2.1) A ft (fe,a:) = (2k + l)hb + h 2 

When k = 0, we get: 
Corollary 2.2. 

.For any 167, £/iere exist C and ho > 0, such that, for any h G]0, ho], there exists 
$g G C C °°(M), $g ^ , such that 

\\ffh$h _ X h (x)^\\ < Ch 17 / 8 \\^\\ , 

where 

X h (x) = hb + h 2 ^-, 

This corollary gives immediately the upperbound in Theorem 11.11 

The proof of Theorem l2.1l is long, so we will split it in different steps in the next 
subsections but let us previously discuss the interpretation of the coefficients. 

2.2. Geometrical interpretation of the coefficients. 

The term 

(2.2) (2k + l)hb + (fc 2 + k) R 

in the right-hand side of (|2.1I) has a natural interpretation. It depends on whether 
R is zero, positive or negative and in all three cases can be described in terms of 
eigenvalues of the associated magnetic Laplacian with constant magnetic field (Lan- 
dau operator) on the corresponding simply connected Riemann surface of constant 
curvature. 

For R = 0, it is a well-known fact that the eigenvalues of the magnetic Laplacian 
with constant magnetic field b on the flat Euclidean plane 

are given by the Landau levels 

X h (k) = (2k + l)hb, keN. 

In the case R is negative, we consider the hyperbolic plane H, which is realized 
as the upper halfplane 

U = {(x,y) eK 2 :y>0}, 
equipped with the Riemannian metric 

dx 2 + dy 2 
■9 = 2 • 

y 



(2k 2 + 2k + l)^+ l -(k 2 + k) R(x) 
Abo 2 v ' w 



4 



BERNARD HELFFER AND YURI A. KORDYUKOV 



Then the scalar curvature is a negative constant R = — 2. The magnetic Laplacian 
with constant magnetic field b on H is given by (see, for instance, pQ) 

This operator first appeared in theory of automorphic forms, where it is known 
as the Maas Laplacian. Its spectrum in i 2 (H) was studied by Elstrodt in [2J. It 
consists of absolutely continuous and discrete parts. The absolutely continuous 
part is given by 

<7 ac (H h ) = {h 2 b 2 + ^,+^{. 

The discrete part is empty if < hb < \ and, if hb > \, consists of a finite number 
of eigenvalues of infinite multliplicity (hyperbolic Landau levels) given by 

X h (k) = (2k + l)hb -h 2 (k 2 + k) , k € N , k < hb - - . 

It is clear that the last formula agrees with (|2.2j) . 

Finally, in the case R is positive, we consider the two-dimensional sphere 

S 2 = {{x, y, z) e R 3 : x 2 + y 2 + z 2 = I}, 

equipped with the Riemannian metric induced by the standard Euclidean metric 
in R 3 . Then the scalar curvature is a positive constant R = 2. In this context the 
magnetic Laplacian is constructed as the Bochner Laplacian, acting on sections of 
a Hermitian line bundle C with a compatible connection V over S 2 . The magnetic 
field is the curvature 2-form B of V. It is constant if B is a scalar multiple of the 
volume 2-form dx g : 

B = s dx g , s € R . 

The construction of a Hermitian line bundle C with a compatible connection 
V such that the associated magnetic field is constant is a particular case of pre- 
quantization in geometric quantization |16j . Such a line bundle exists if and only 
if s = n/2 for some n £ Z. For s = ra/2, the corresponding line bundle (£„,V„) 
can be described as a complex line bundle associated with an S^-principal bun- 
dle S 3 —> S 2 called the Hopf fibration and the character \ n : S 1 — > S 1 given by 
Xn(u) = u n ,u € S . In physics literature, (£„, V„) is a well-known Wu-Yang mag- 
netic monopole [26] . which provides a natural topological interpretation of Dirac's 
monopole of magnetic charge g — nh/2e. 

The magnetic Laplacian H n acting on sections of C n is defined as 

H n = v;v„ . 

Its spectrum was computed in [27] (see also [13 HI]). It consists of a countable set 
of eigenvalues (spherical Landau levels) given by 

(2.3) ^\n\(2k + 1) + k 2 + k , keN, 

with multiplicity \n\ +2fc+l. The corresponding eigenvalues are known as monopole 
harmonics. The formula (|2 . 3[) agrees with (|2 .2[) if we take h = 1/n and 

h h = ^v;v„ . 

n z 

Finally, let us mention the paper [3] , which states that the three types of magnetic 
Laplacians described above are integrable in some sense and provides the complete 
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description of their spectra in the same way as it was done by Schrodinger for the 
harmonic oscillator. 

2.3. Expanding operators in fractional powers of h. 

The approximate eigenfunctions € C£°(M), which we are going to construct, 
will be supported in a small neighborhood of the point x S 7 appearing in the 
theorem. We will consider some special local coordinate system with coordinates 
(s, t) in a neighborhood of 7 such that 7 corresponds to t = 0. We will only apply our 
operator on functions which are products of cut-off functions with functions of the 
form of linear combinations of terms like h v w(s, h^ x / 2 t) with w in C°°{S 1 ) <g)<S(Rt). 
These functions are consequently 0(h°°) outside a fixed neighborhood of 7. We will 
start by doing the computations formally in the sense that everything is determined 
modulo 0(h°°), and any smooth function will be replaced by its Taylor's expansion 
at 7. It is then easy to construct non formal approximate eigenfunctions. 

Choose a normal coordinate system in a tubular neighborhood U of 7 with 
coordinates X = (Xq,Xi) with Xq = s and X\ = t. Here s e [-L/2, L/2) = 
S\ = K/iZ is the natural parameter along 7 (L is the length of 7), 7 is given by 
the equation t = 0, and t G (— £o,£o) is the natural parameter along the geodesic, 
passing through the point on 7 with the coordinates (s, 0) orthogonal to 7. 

It is well known that in such coordinates the metric g has the form 

g = a(s, t)ds 2 + dt 2 , 

where 

a(s,0) = I. 

So we can write 

(2.4) g o(s,t) = l + a 1 (s)t + a 2 (s)t 2 +(D(t 3 ), g 01 (s,t)=0, g n (s,t) = l. 
In particular, we have for the first coefficient of the inverse matrix {g 13 ) 

(2.5) g m ( S ,t) = I - ai (s)t - (a 2 (s) - ai ( S ) 2 )t 2 + 0(t 3 ). 
It is known (see (|A.1[) and (|A.2|) in the appendix) that 

1 2 

<2l = — 2k, d2 = — —it + K , 

where k is the mean curvature of 7. 
Let us define 

\g\ = det{ 5ij } = gooSn - 3oi > 
A = A ds + Aidt, 

and 

dA 1 dA 



(2.6) B = b{s,t)y/\g(s,t)\dsAdt, b(s,t)y/\g{s,t)\ 
The external unit normal vector to 7 has the form 

and its natural extension N in the neighborhood of 7 is 

- d 



ds dt 
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Without loss of generality, we can assume that b(x) > for any x £ M. Thus, we 
have 

db d 2 b 

(2.7) 6(«,0) = 6o, ^(s,0)=0, /3 2 (s) = ^(s,0)>0. 

We can assume that (after a gauge transformation) 

Ax(s,t) = 0. 
Using (EU) and ([221), we obtain that 

(2.8) A (s,t) 

= ai(s) - b t - ^ ai {s)b t 2 - i (#,(«) + b (a 2 (s) - ^ai(s) 2 )^ t 3 

+ 0(t*), 

as t — > . Our constructions will be local, in a neighborhood of i £ 7. So, using 
gauge invariance, we can take a± = 0. 
We have 

H h = -=L= E V " ( v / R^)l5 Q/3 (^)V^ , 

VlSWI 0<«,/5<l 



where 



V£=i/i^+A*P0, a = 0,1 



or, taking into account (|2.4|) . 

rv'^M 2 _ t> : ,,,, , 

at 2 at 



(2.9) if' 1 = .g 00 ( S ,t)(V(}) 2 - + t7»T°V§ - h 2 T^ 



where 

r° = - 1 V -|- f ^M/ a (^)) , a = 0,1 

By and ((2~5|) . it follows that 

t^O 1 ^ ( /TZ7v^\ J)0 r v\\ 1/ 



1 a 



(/Rx^ 00 (X)) = --oi(*)* + G(i 2 ) , 

( v 7 ^) = ^oi(«) + (02(a) - ^«i00 2 ) t + o{t 2 ) . 



We now move the operator H h into the Hilbert space L 2 (U, ds dt), considering 
the operator 

(2.10) H h -^IgWfHhlgiX)]- 1 '* . 
By (EH), we get 

(2.11) H h = g Q0 (s, t)(V%) 2 + (V^ 1 ) 2 + ihT°V% + ihT 1 ^ , 
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where 

Vg:= \g(X)\^V h \g(X)\-^ 
(2 " 12) Vl:=ih\ 9 (X)\^\g(X)\-^ 

= ih \j t - - - \*i{sf)t + o(t 2 ) 

Now we use the scaling t = h}/ 2 ti and expand the operator H h in power of h 1 / 2 . 
By straightforward computations, we obtain that 



H h = hH 



/!.() 



where 
with 



H h,o = Pq + h 1 ' 2 P 1 + hP 2 + h 3 ^ 2 U 3 {h) . 





d 2 

dt 2 4 


-&o*i> 


Pi = 


- 26 *i 


(4 




d 2 


3i 



d 1 

+ \ (a2(s) - ^ai(s) 2 

and TZs(h) a second order differential operator whose coefficients are formal power 

series of the form S^lo h-*' 2 aj{s,t\). 

2.4. Construction of approximate eigenfunctions. 

Let us fix x € 7 and k e N. Without loss of generality, we can assume that 
x corresponds to <E S\. We introduce a self-adjoint second order differential 
operator Q h '° in L 2 {S\ x M), which is close to H h '° near {0} x R, and construct 
approximate eigenfunctions tp h of the operator Q h, ° in the form 

fh = <Po + h 1/2 ipi + h(p 2 , 

such that 

\\(Q h '° - \(h))vh\\ =o(h 3 / 2 ), 

with X(h) of the form \(h) = Ao + /i 1/,2 Ai + /iA2 . (Observe that such approximate 
eigenfunctions don't exist for the operator H h '°.) The fact that the operator H h '° 
is close to Q h '° near will allow us to construct approximate eigenfunctions for 
H '°, using appropriate /i-dependent cut-off functions, and complete the proof of 
Theorem [2H1 

The operator Q h '° is defined by the formula 

Q h '° := L + h 1/2 Lx + hL 2 , 
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where 



Lo = 


d 2 

ml 4 




Li = 


- 2Mi 




L 2 = 


d 2 

ds 2 


3? 8 1 
-foi(«)^6o*? + 560^2(0) 



- \a 2 (0)blt{ + l - ai {0) 2 blt{ + I (a 2 (0) - Ja^O) 2 

+ j a'i(*)M? + ^ai(s) 2 ^i " ^«i( s ) 2 ■ 

Now we write formal expansions in powers of h 1 ^ 2 : 

00 00 

3=0 j=o 

and express the cancellation of the coefficients of hp' 2 in the formal expansion for 
(Q h <° - X(h))cp h for .7 =0,1,2. 
Step 1 

At the first step we get 

(2.13) L tpo = \ (p . 

We have 

Sp(L ) = {nm = (2m + 1)6 : m e N} . 
The eigenfunction of Lq associated to the eigenvalue fi m is 

i> m (h) = n-V%/ 2 H m (bX%)e- b °^ 2 , 
where H m is the Hcrmite polynomial: 

H m (x)=(-l)V^(e- 2 ). 

The norm of ip m in L 2 (M.,dx) equals the norm of H m in L 2 (E, e~ x dx) , which is 
given by 



\\H m \\ = ^/2™m!V^. 
For the given fixed integer k > 0, we take a solution of (|2.13l) in the form 

A = (2k + l)b , Po(Mi) = Xo(s)i/Jk(h) ■ 
We recall that the Hermite polynomials satisfy 

H m+ i(x) = 2xH m (x) - 2mH m _ 1 (x) , 

and 

H' m (x) = 2xH m {x) - H m+1 . 

We also have 

2xH m =H m+ i + 2mH m _i, 

Ax 2 H m =H m+2 + (4m + 2)H m + 4m(m - l)H m - 2 , 

8x 3 H m =H m+ 3 + (6m + 6)H m+ i + 12m 2 H m _i + 8m(m - l)(m - 2)fl m _ 3 . 
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l6x i H m (x) =H m+i (x) + (8m + 12)H m+2 (x) + 12(2m 2 + 2m + l)H m (x) 
+ 16(2m 2 - 3m + l)mH m - 2 (x) 
+ \Qm(m — l)(m — 2)(m — i)H m -^{x) . 

Using these identities and the orthogonality of the Hermite polynomials, we get 

1 , , , \ 1 1 , , , \ k + 1 

{hipk-i , ipk) = — , n TT U7 {tiipk+i,ipk) 



W — — 26 l/2' ||if fc ||^-— ^ 6 l/2 > 

1 11 3 



— (tfyk+3,V%> = 4?( fc + 3 )( fc + 2 )( /:: + 1 )'^n^( i i ? /'fc+i^fc) = ^(fc + i) 2 



We will also use several times the following identities: 

' ;((4D t2 -liti)i> k ,ip k )=0, 



\\H k \\ 2 \ z 12 2 

' ;(l4^ k+3 + (3k + 3)t|Vfe+l - efc 2 ^^.! - -k(k - l)(k - 2)i^*-3, Vfc) 



(l5fc 2 + 15fc+^ 
&o /2 V 2 



and 



^2 (*2 + (3k + 3)^ fc+ i - 6fc 2 V> fc -i - - l)(fc - 2)^ fe _ 3 ) , V'fc) 

(2fc + l) 



\Hi„ 

3 



1/2 



Step 2 

At the second step, we obtain 

(2.14) (L - \o)(fi = Xipa - Litpa , 

or, more explicitly, 

-||+& 2 < 2 -(2fc + l)6o) <p! 



=AiXo(s)Vfe(*i) + «Xo( s ) & o /2 (' i /'fc+i(^i) + 2fc^fc-i(ii)) 



+ ai(s)xo(s)&o /2 (^?Afe+3(ti) + f|*+|j ipk+i(ti) 

+ jfc 2 Vfc-i(*i) + ^fc(fc - l)(fc - 2)Vfc-a(fi)) • 
The equation (|2.14[) has a solution only when 
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and this solution can be taken in the form 

i 

ViOMi) =— 575X0 {s)[ipk+i{ti) - 2kip k -i{ti)] 
2b 



+ a 1 (s)xo (5)^72(^^+3^1)+ |j ^fc+i(*i) 

- |fcVfc-i(ti) - ifc(fc - l)(fc - 2)^ fc _ 3 (ti)) • 

Step 3 

At the third step, we have 

(2.15) (L - X )ip 2 = A 2 </?o - L 2 ip - Litpt . 

This equation has a solution if and only if, for any s, its right-hand side is orthogonal 
to ipk(h)- 

First, let us compute L 2 (^o : 

L2V0 = - Xo(s)Mh) + yai(*)Xo(*)Mi^fc(*i) + ^2(0)xo(s)M^(*i) 

- ^ 2 (0) X o(s)6o*iV'fc(ii) + ^i(0) 2 Xo(s)^i<M*i) 

+ I ( a 2 (0) - 7ai(0) 2 ) Xo(*)^fc(ti) + ^ai(s)xoW&ot 2 ^(ti) 



2 V 4 
+ ^ai(s) 2 Xo(s)&^^ fe (ii) - -^i(s) 2 XoO#fc(*i) ■ 
Multiplying by V'fc and integrating with respect to t\ gives: 

\\H k \\2 y L 2 (po(s,t 1 )i>k(h)dt 1 

= x'o(s) + ai(*M*)i (f* + J) + ^(2fc 2 + 2fc + l)ft(0) X o(«) 

- (fc 2 + fc)a 2 (0) X o(s) + J(fc 2 + k) ai (0) 2 X o(s) 

+ j(2k + l)oi(*)xo(«) + {^k 2 + ai ^ Xo ^ ■ 

Next, let us compute L\ip\\ 

Li<pi ^bl /2 Xo(s)[tiipk+i(ti) - 2ktiip k -i{ti)] 

1 ( 3 3 \ 

- a'i( s )xo(s)^ /2 i(^iiV'fe+3(ii) + ( g fc + g ) 

- Jfc 2 ^fc-i(ii) - ^(fc - l)(fc - 2)iiVfc- 3 (ii)) 



- a i( s )Xo( s ) & o /2 i(^*iV'fe+3(ii) + (Jfc + I J *i^fe+i(*i) 

- |fc 2 tiVfc-i(*i) - ^(fc - 1)(* - 2)ti^ fc _ 3 (*i)) 

- ^i(s)Xo(s)&o /2 [*?V'fe+i(ti) - 2fct?Vfc-i(*i)] 
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- ^jfc 2 ^fc-i(*i) - ifc(fc - l)(fc - 2)i?Vfc-s(ti)) ■ 
Multiplying by ipk and integrating with respect to t\ gives: 
1 



1 



=Xo( s ) ~ a'i(s)xo(s)i (^k + |j - ai(s)xo(s)i Qfc + 



Thus, multiplying the right-hand side of (|2.15p by ipk an d integrating with respect 
to ti we obtain 

1 



\H k V z 



=A2Xo(s) - ,, „ ll2 / L 2 ipo(s,ti)ipk{ti)dt 1 - ^ / Li<£i(s,ti)^ fc (ii) d*i 



it 



=A 2X0 ( S ) - + 2fc + l) X o(s) + (fc 2 + fc) L(0) - Jai(0) 2 ) *„(«) . 

So the solvability condition for ()2.15p holds if we put 

A 2 = ^^(2A; 2 + 2k+ 1) - (fc 2 + fc) U 2 (0) - Jai(0) 2 ^) , 

and, in this case, there exists a solution y> 2 of (|2.15p . 

Thus, we arrive at the following proposition. 
Proposition 2.3. 

For any xo € C£°(K)> there exists an approximate eigenfunction ip^ of the operator 
Q h, ° in the form 

(2.16) ip h = ip Q + h 1/2 tpi + hip 2 , 
where tpo *s given by 

(2.17) Ms,h) = Xo(s)Mh) , 
ipi and ip 2 have the form 

(2.18) tp^s, h) = A {s, h) Xo (s) + At(s, ti)xo(«) , 

(2.19) y 2 {s, tt) = B (s, h)xo(s) + ti)Xo(*) + B 2(s, *i)Xo(*) . 
where Aq, A\, Bq, B± and B 2 are some smooth functions, such that 

\\(Q h < -\(h))tp h \\ = O(hV 2 ), 

with 

X(h) = A + h\ 2 , 

where 

A = (2k + l)bo , A 2 = ^^(2fc 2 + 2fc + 1) - (fc 2 + fc) (a 2 (0) - Jai(O) 2 ^) . 
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Observe that X(h) — h 1 \ h (k,x), where X h (k,x) is given by (|2~T 



Now we consider the function as constructed in (|2.16[) in the previous propo- 
sition but with xoj depending on h: 

X0 (s) ■= E h (s) = ch-^ 2 X ( S )exp (-^3) ■ 

Here f3 e]0, l/2[ is some constant, which we will choose later, and \ € C£°(M) is a 
cut-off function, which equals 1 in a neighborhood of 0. 
Observe that 

(2.20) \\s m E h (s)\\ = (V^ S 2m exp (-^) ds\ ' = C^ m , 
and, furthermore, 

(2.21) \\s m E' h {s)\\ = C 2 h^ m - X \ \\s m E'l{s)\\ = C 3 ^ (m ~ 2) . 
We have 

(Q h '° - X(h))<f> h = h^ 2 {L 2 ^ + L W2 - \ 2ipi ) + h 2 (L 2 ^ 2 - X 2 ip 2 ) . 

Since L\ is a first order differential operator in s and L 2 is a second order differential 
operator in s, by (|2.18p and (|2.19l) . it follows that 

L 2 y x + L x <p 2 = F E h + F 1 E' h + F 2 El + F Z E'£ , 

where Fj, j = 0, 1, 2, 3, are some functions. Using (12.20)) and (|2.2ip . we obtain the 
existence of constants C > and ho > such that 

||(<2 M - \(h))<5>h\\ < Ch^ 2 -^\\^ h \\ , Vh e}0,ho] ■ 

Then we have 



H h '° - Q 
h 



h,0 



\b t\{p 2 {s) - fc(0)) - l(a 2 (s) - a 2 (0))b 2 4 + Ua.is) 2 - ai {Q) 2 )b 2 t\ 



1 



a 2 (s)--ai(sY - U 2 (0)-- ai (0) 



1 



which immediately implies that 





\\(H h '° - 


Q hfi )®h\ 


< C7i 7o(/3) ||$/> 


with 










iM 


= min(l - 


- /3, 3/2 - 2/9) , 


and 










m ,,o_ 


■ \(h))* H 


< CKiW\\* h 


with 










7C8) = 


= min(l -f 


£,3/2-3/3). 


Setting 









1 



we obtain the existence of a constant C such that 

\\(H h <° - \(h))$ h \\ <C0 B \\$ h \ 



h 3 ^ 2 1l 3 (h) 
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3. Lower bounds 

In this section, we will prove the lower bound for the groundstate energy Xo(H h ) 
of the operator H h . As above, we write 

PLq := inf /3 2 (s) . 

Theorem 3.1. 

There exist C and ho > 0, such that for any h €]0, ho] 

\o(H h )>hb + h 2 ^-Ch 19 / s . 

First, we recall a general lower bound due to Montgomery [22]. Suppose that U 
is a domain in M. Then, for any u £ C%°(U), the following estimate holds: 



(3.1) \\(ihal + A)u\\(j> h 



b\u\ 2 dx & 



This fact is an immediate consequence of a Weitzenbock-Bochner type identity. 

From (|3.1[) , it follows that we can restrict our considerations to any sufficiently 
small neighborhood of 7. We denote by the Dirichlet realization of the 
operator H h in L 2 (£l, dx g ) . 

The estimate (|3.ip implies that 

ThH% + (1 - rh)hb < OKTKh- 1 . 
Taking r = h~ x / 2 , we obtain 

h 1,2 {H h D -hb + h 1/2 b) <H h D , 0</i<l. 
Consider the Dirichlet realization P^ of the operator 

(3.2) P h = H h -hb + h^ 2 (b-b ), 
in L 2 (fl, dxg). Then we have 

(3.3) hbo + h l / 2 Xo{P^)<Xo{H h D ). 

Therefore, the desired lower bound for Ao(-ffjj) is an immediate consequence of the 
following proposition. 

Proposition 3.2. 

There exist Co and ho > such that 

* (Pb)>h a '*£--C h 1B '\ he]0,h }. 



To prove Proposition I3.2[ we will follow the lines of the proof of [12J Theorem 
7.4]. First, we observe that the upper bound in Theorem 11.11 and (|3.3p imply an 
upper bound for Xo(P^): 

(3.4) Xo(Pri) < /i 3/2 + Ch 13 / S . 

46 

Denote by ut an eigenfunction of associated with Xo(P^): 

(3.5) P*u h = X (Pk)u h . 

By a straightforward repetition of the arguments of [12], we can easily show the 
following analogue of Lemmas 7.10 and 7.11 in |12pl 



^There are a few inaccuracies in |12| . concerning Lemma 7.11. For the erratum, see 
http: / / www.math.u-psud.fr / ~helffer /erratuml64II.pdf 
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Lemma 3.3. 

For any real k > 0, we have 



\\d(x r/ ) k u h \\ LHn . dXg) < c kJ (h k / 2 + h^ k+i y s )\\u h \\ LHn4Xg) . 

For any a = 0, 1 and any k > 0, we have 

As above, take normal local coordinates X = (X , X±) = (s, t) near 7 such that 
7 corresponds to t = and assume that they are defined on Q. Thus, we have 

b(X) = b Q + ^ 2 ( s )t 2 +0(t 3 ). 

and, for the metric coefficients, 

g 00 (s,t) = l + a 1 (s)t + a 2 (s)t 2 +0(t 3 ), g ol (s,t) = 0, g n (s,t) = l. 

We can choose a magnetic potential A such that 

1 



A (s,t) = ai(s) - b t - ^ai(s)M 2 + 0(t 4 ), A X (X) = . 



We have 



p h D = ]T - / ^v h a Vbm\g af) (x)v h ,-hb(x) + h^(b(x)-b ), 

o< a ,0<i V IflWI 
so its quadratic form is given by 

{P%u,u)=f 9 aP (X)V h a u(X)V^u(X)^giX)dX 

n 0<a,/3<l 

-hi b{X)\u{X)\ 2 JgJX)dX + h x l 2 [ (b(X)-b )\u(X)\ 2 ^gJX)dX. 
Jn Jn 



Note that 



Pd > 0. 



Now we move the operator into the Hilbert space L 2 (f2,dX), using the 
unitary change of variables v = \g(X)\ 1 ^ 4 u. For the corresponding operator 

P^\g(X)\^P^\g(X)\-^ 

in L 2 (tt, dX) , we obtain 

(Pkv,v)=[ J2 9 a ^X)(v h a -Uh\g{X)\-^\g(X)\)v{X)x 

-/j f b(X)\v(X)\ 2 dX + h}' 2 f (b(X) - b Q )\v(X)\ 2 dX . 
Jn Jn 



Put 

q ^ = l E KvS-^lffWr^lffWl)^)! dX-hb Jjv(X)\ 2 dX. 



Then we have 

(3.6) \{P*v,v) - q(v) 



SCHRODINGER OPERATORS WITH MAGNETIC WELLS 



15 



< 



^-lih\g(X)\-^ s \g(X)\jv(X) 



dX 



+ C\h / t 2 \v{X)\ 2 dX + C 2 h 1 l 2 / t 2 \v{X)\ 2 dX. 



Consider the Dirichlet realization P^ od D of the operator 

PLd = (ih§~ s - hot) - h 2 ^ - hbo + h 1 ' 2 Q/3 2 ( S )t ; 
in the space L 2 (£l, dX). So its quadratic form is given by 



(P* od v,v) = q™ d (v) + hV 2 Q/? 2 (^ 2 ) HX)\ 2 dX, 

dX-hbo [ \v(X)\ 2 dX. 



ih— - b t ) v(X) 



dX+h 2 



where 
q mod {v) = 
So we have 

(3.7) \{Pi,V,v) - (PL d V,v)\ < \q(v)-q mod (v)\ 

+ Y1 I t 2 \V h a v{X)\ 2 dX + h [ t 2 \v{X)\ 2 dX + h 1 ' 2 [ \t\ 3 \v(X)\ 2 dX 
a Jn Jo, Jn 

Finally, we have 

(3.8) \q(v)-q mod (v)\ 

<C(q(v)) 1/2 h( [ t 2 \v(X)\ 2 dx) +( [ t 6 \v{X)\ 2 dX 



1/2 



By Lemma l3~3l for any real k > 0, there exists C k > such that 

(3.9) \\\t\ k u h \\ LWxg) < C k (h k / 2 + h^ k+i y s )\\u h \\ L2 ^ dXgh 
and, for any a = 1,2, 

(3.10) \\\t\ k V h a u h \\ LH n, dXg) <C fe (^ fe+1 )/ 2 + M 3 ' £+5 )/ 8 )||^|| i2( a (i:C3 ). 

Put v h = {giX^^Uh. By and (|3~TU|) . for any real k > 0, there exists 

Ck > such that 

(3.11) \\\t\ k v h \\ LH ^ dX) < C k {h k l 2 + h^ k+i y & )\\v h \\ L ^ dX)l 
and, for any a = 1,2, 

(3.12) 



l 2k 



d 

ih— - b t ) Vh(X) 



2 \ i/a 
dX 



\t\ 



2 k 



h^iX) 



2 \ 1/2 

dX 



<c k (h( k +^ + h^ +5 y s )\\vd W xy 

The estimates (|3.11[) and (|3.12[) allow us to show, first, using (|3.5[) . f|3 .4[) and 
HI), that, 

< Ch*' 2 \\ 

v h\\ 2 L2(n,dX)i 
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then, using p.8[) . that 

\q{v h )-q mod {v h )\<Ch 2 \\v h f L2{w , 
and finally, using (|3.7I) . that 

(3.13) (PLdV h ,v h ) < (A (P^) + C/ i 15 / 8 )||^||l 2(aidJf) . 

Let x be a function from C£°(f2) such that \ = x(t) an< i X = 1 in a neighborhood 
of 7. By (|3.11|) and (|3. 12|) . it follows that, for any k E N there exists Ck > such 
that 

<9y , 

(3.14) ||(1 - x)vh\\mn,dx) + \\-Q^Vh\\L^(n,dX) < C k h '\\v h \\ L 2 (n4X) . 

and 
(3.15) 



(i-x(*)) 



z/i— - 6 i ) Ufc(X) 



2 \ Va 



(i-x(*)) 



1/2 



Using (|3.14[) and ()3. 15|) . it is easy to check that, for any k > 0, there exists 
Ck > such that 

(3.16) \{PLdiXVh),XVh) ~ (PL*>h,v h )\ < C k h k \\v h f. 

Consider the self-adjoint realization of the operator P^ od in L 2 (M. 2 , dX). We will 
keep the same notation P^ od for this operator. Put Wh = x\9{^)\ X ^ u h- By (|3.13[) 
and (|3.16l) . it follows that 

(P* od w h ,w h ) < (A (P£) + Co/i 1578 )!! 
This immediately implies that 

(3.17) X (P^ od )<X (P^) + C h 15 ^. 
Consider the operator 

Recall that the eigenvalues of the Schrodinger operator with constant magnetic field 
and positive quadratic potential in 1" can be computed explicitly. More precisely 
(see for instance [23 Theorem 2.2]), the eigenvalues of the operator 

2 / cv 1 \ 2 

+ K n x 2 + 2K l2 xy + K 22 y 2 



Hb,K 

are given by 
where 



4 " \ hV 



1 % + \ hX 



X ni n 2 = (2ni + l)si + (2n 2 + l)s 2 , ni,n 2 £ N. 



si 



S2 



1 

1 

71 



^ + 6 2 -[(^ + & 2 ) 2 -4^] 1/2 



K 



*,2\2 



4d*] 1/2 



1/2 



1/2 



and 



tjc = Tr if = Jf u + K 22 , = dot K = K n K 22 - K 2 2 . 
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Applying this formula to the operator Pq 1 , we obtain that its eigenvalues have the 
form (t K = |/io, d K = 0): 

A ni n 2 = (2ni + l)si + (2ra 2 + l)s 2 - hb , ni,n 2 eN, 

where 

si = 0, 

and 

V 2 , , i x 



J 2o 



s 2 = h 

For the lowest eigenvalue Aq(Pq) of Pq , we obtain 
(3.18) X (P^ = ^h 3 / 2 +O(h 2 ). 

Observing that A (P^ od ) > A (P ' 1 ) and comb ining (|3 . 1 T[) and (|3.18[) . we immedi- 
ately complete the proof of Proposition [ 



4. MlNIWELLS 

4.1. Main statement. 

Fix some k £ N and additionally assume that there exists a unique minimum xq G 7 
of the function 

(4.1) Vfc(a:) := (2k 2 + 2k + 1)^M + I ( fc a + fc ) R{x) 

on 7, which is nondegenerate, that is satisfying, for all x € 7 in some neighborhood 
of xq, 

(4.2) Cd(x,x ) 2 < V k (x) - V k (x ) < C- l d{x,x f . 

The purpose of this section is to give the following more precise construction of 
approximate eigenvalues of the operator H h . 

Theorem 4.1. Under current assumptions, for any natural j , there exist Uj k € 
C%°(M), C jk > and h jk > such that 

= S hh + °h,h,k(h) 

and, for any h G]0, hjk], 

\\H h u%- ^u%\\<C ]k h^\\u%\\, 

where 

(4-3) ^jfc = Mj,fc,o/i + f-i 3 ,kAh 2 + Hj, kfi h b/2 , 

with 

fij,k,o = (2k + 1)6 , Hj,kA = Vk(x ), 

and 

Mi.M = -^14"(xo) 1 / 2 /3 2 ( 2 ;o) 1 / 2 (2fc + l)^ 2 (2j + 1). 

Here and below the derivative means the derivative with respect to the natural 
parameter on 7. 
As above, denote 



(J, := inf /3 2 (x). 



When k = 0, we get: 
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Corollary 4.2. 

Assume that there exists a unique minimum point xq G 7 of the function fi% on 7, 
which is nondegenerate : 

H 2 := %(x ) > 0. 

For any natural j, there exist € C^°(M), Cj > and hj > such that 
and, for any h e]0, hj], 

\\H h u$-$v.$\\<C j h 11 ' i \\u% 

where 

Vj=hb + h — + h 1 (2j + 1). 

4& 46^ 

Corollary 4.3. 

Under the assumptions of Corollary \4-l\ for any natural j, there exist Cj and 
hj > ; such that for any h g]0, hj] 

Xj(H h ) <hb + h 2 ^ + ^ 2 { ^fj 2 (2j + 1) + Cjh n ^ . 
4b o 46g 7 

Remark 4.4. We conjecture that under the assumptions of Corollary 14. II 

Xo(H h ) = hbo + h ^ + h */* ^1^1 + o(h^) . 
46 

4.2. Expanding operators in fractional powers of h. 

The approximate eigenfunctions £ C°°(M), which we are going to construct, 
will be supported in a small neighborhood of 7. As in Section l2~3l we will use the 
normal coordinate system (s,t) in a tubular neighborhood U of 7 with coordinates 
X = (Xq,X\) with X a — s and X\ — t. We make use of notation of Section [231 
and need to consider further terms in the asymptotic expansions of that section. 
Thus, instead of (|2.4[) . we write 

g Q0 {s,t) = l + a 1 {s)t + a 2 { S )t 2 +a 3 (s)t 3 +(D(t 4 ), 

(4.4) 

gu(s,t) = l, g 01 (s,t)=0. 



Then (|275|) takes the form 

(4.5) g 00 (s,t) = 1 - ai (s)t - (02(a) - ai(s) 2 )i 2 

- (a 3 (s) - 2ai(s)a 2 (s) + ai(s) 3 )i 3 + C(t 4 ) , * -> . 
For the components of the magnetic potential, we assume that 

Ai(s,t) = 0, 

and, instead of (|2.8j) . we get 

(4.6) A (s,i) ="1 -M- jai(s)M 2 

1 

~ 6 



182(a) + b (a 2 (s) - - ai ( S Y) t A + 0(t 4 ) , t 



where ai is some constant. Without loss of generality, we can locally (after a gauge 
transformation) take a\ = 0. 
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Wc have 



r° = - iaiOO* + (-\a' 2 (s) + a^a'^s)) f + 0(t 3 ), 
(4.7) r l =~ ai {s) + ( fl2 ( s )-i ai ( s ) 2 ^ 

+ (|oa(«) - |ai(*)02(«) + ^i(*) 3 ) ^ + 0(t 3 ) . 



We move the operators into the Hilbert space L 2 (R") equipped with the Eu- 
clidean inner product. For the operators Vq and defined in (|2.12l) . we obtain 
the following expansions: 



Vg =Vg - ^ha' 1 (s)t--h(a' 2 (s) - ai {a)a[{a))^ + 0(t 3 ) 
(4.8) = - /iA -oi(s) - -(a 2 (s) - ^ai(s) 2 )i 



We will assume that the minimum point Xq £ 7 corresponds to s = 0. Then the 
condition (|4.2j) implies that 



®W + 2fc + 1) - (fc 2 + fc) (a 2 '(0) - ia!(0)oi'(0) - ^(O)) 2 ) > . 



We also have 



(4.9) + 2fc + 1) - (fc 2 + fc) (a 2 (0) - iai(0)oi(0)) 



Now we use the scaling t = h x / 2 ti, s — /i 1 / 4 s 1 . Similarly to Section |2~51 we will 
only apply our operator on functions which are a product of cut-off functions with 
functions of the form of linear combinations of terms like 1 / 4 s, h~ 1 / 2 t) with 

w in C 00 ^ 1 ) ®S(R t ) supported in a small neighborhood of (0, 0) £ 5 1 x R t . These 
functions are consequently 0(h°°) outside a fixed neighborhood of (0,0). We will 
start by doing the computations formally in the sense that everything is determined 
modulo 0(h°°), and any smooth function will be replaced by its Taylor's expansion 
at (0,0). 
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(4.10) 



From (|4.5p . (|4.6|) . (|4.7|) . ()4.8j) . we derive the following expansions: 
10) 

:i + itfl*-?- - iai(0)6oft*2 - 7a' 1 (0)6 / l 5/4 Sl ^ 
asi 4 4 



^ 3 / 2 



f/3 a (0) + 6o(a 2 (0) - Jai(0) 2 )) t? - -4(0)*! - g4'(0)Mf* 



2 + 2 



/3 2 (0)+6 (a 2 (0) - ioi(0)ai(0) 



sit? 



4 
^ 2 



^a' 1 '(0) Sl t 1 --af(0)&o^ : j 

1 fe(0) + 6oK(0) - i ai (0)a'/(0) - ^(4(0)) 2 )) s 2 t 



2 + 3 

1 



\<W\t x - ^a[ v (0)b o sitt - i (a 2 (0) - ai (0K(0)) * 



- 0{h 9 ' 4 ) . 



(4.11) 



h 1 ' 4 



fi __ /iai(0 )--^/v 1 (o) Sl 

i^O)--^) 2 )^^)^ 



ii(a 2 (0) - ai(0K(0)) Sl ii + ^<(0)s? 



|(5a3(Q)-5a 1 (0)a2(0) + ia 1 (0) 3 )/7 



i (a 2 '(0) - Ol (0)oi'(0) - K(0)) 2 )^i + ^4 v (o)4 



0(h 9/4 ). 



(4.12) 



<? 00 =1 - axityh^hx - a[(0)sxh 3 /Hx 



h 

h 5 ' 4 



-(a 2 (0) - ax(0) 2 )tl - ~a'{{Q)s\tx 



--a'"(0)s\tx - (a' 2 (0) - 2a 1 {Q)a' 1 {0))s 1 t\ 



-(a 3 (0)-2a 1 (0)a 2 (0) + a 1 (0) 3 )i? 



~(4'(0) - 2 Ol (0K(0) - 24(0) 2 )s 2 i 2 - ±a{ v (0)s$tx] + O^ 4 ) 



(4.13) 



=/ l 1 / 2 



-4(o)*i 



ft 3 / 4 



jai(0)«iti 



-4(0) + ai(0)oi(0)U?-4oi"(0)^i 



0(/i 5 / 4 ) . 
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r 1 = 



I ai(0 ) + h^a'Msr + h 1 / 2 [L(0) - i ai (0) 2 ) t x + ±a?(0) fl ? 



ft 3 / 4 



(4.14) 



12 



«i"(0)s? + (4(0) - ^(0)4(0)) sxh 



-a 3 (0)--a 1 (0)a 2 (0) + -a 1 (0) 3 W? 



+ \ (4 (0) - ^(0)4(0) - (4(0)) 2 ) sjh + ^a[ v (p)st] + 0(h 5 ^) . 

Next, we make the partial Fourier transform F Sl ^ HJl in the Si-variable and the 
translation 

t 2 = t 1 + h 1/i< ^-, a 2 = ai, 
oo 

and expand the operator H h defined by (12.10)) in powers of h x l A . Using (12~TD) . 
(|4TT0l) . (|4TTTj) . (|4TT2"1) . (|4TT3|) . ((4TT4|) . after routine computations, we obtain 

H h = hH hfi , 

where 

H h >° = Po + h l ' i P 1 + h l ' 2 P 2 + h 3 / i P 3 + hP 4 + /i 5 / 4 P 5 + h a / 2 P 6 + 0(h 7 ^) , 
with 

9 2 2 2 
= - 7772 + bot 2 , 



P 1 =0. 



P2=- g «1 (0)^0*2. 

P3=^ai(0)6^|D ff2 , 



A =i U(0) - 2b L(0) - |«i(0) 2 )) 6o4 + 5 (oa(0) - ^(O) 2 ) 



+ -ai(0)t 2 4 + ~4(0)&o ^A 2 - 77**2 



/32(0)-6o(a 2 (0)-- ai (0) 2 ) \t\o 2 



:b 



93 

^(0)-26o(4(0)-— Ol (0)oi(0) 



t 4 n 



1 

12 



ja'/(0)6 (2t|A 3 Ar 2 - 3ii|A, 2 ) - ~4(0)t 2 a 2 (eT 2 .D CT2 - 2t 2 A 2 ) 

1 

2 



4(0) - - Ol (0)oi(0) ) 



The formula for Pg is quite long, but we will only need the even part of Pg given 

by 



P+ 



CT) 1 

6b 3 



4(0) -|ax(0)ai'(0)-g(aU0)) 2 



°0 Z 2 
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+ Qo£(0) - ^ ai (0)a'{(0) ^K(O)) 2 ) 



-<(0)6 t 6 2 D t2 - -it: 



Dt 



~ l«i(0) 2 ^ - ^ai(0)(2t 2 A 2 - i) 
4.3. Construction of approximate eigenfunctions. 

Now we look for an approximate eigenfunction ip h of H h '° and the corresponding 
approximate eigenvalue X h as formal expansions in powers of ft, 1 / 4 : 

oo oo 
V H (t2 , (72) ~ J2 ^'^l (*2, <72) , X h ~Y, ■ 

e=o e=o 

In the next steps we will express the cancellation of the coefficients of h e / 4 in the 
formal expansion of (H hfi - \ h )ip h for £ = 0, 1, . . . , 6 . 
Step 1 

At the first step we get 

P (p = Aq^o , 

so we take 

A = (2k + l)b , <Po(t2,<72) = Xo{<72)ipk{h) , 

with the given k. 
Step 2 

The second cancellation gives: 



so we take 



(P - \))<pi = Ai</?o - Pifo = Ai^o • 
Ai=0, <p!=0. 



Step 3 

The next equation reads: 

(Pa - X )if2 = A 2 ^o - P2P0 - P1V1 = -Wo - P2P0 ■ 
The computation of P2P0 gives: 

= - ^a 1 (0)^ /2 X oM(V'fc+3(*2) + (6fc + 6)^ fc+ i(t 2 ) 
+ 12fc 2 Vfc-i(*2) + 8fc(fc - l)(fc - 2)Vfe- 3 (t2)) • 



Thus, we have 
and 



A 2 = 0, 



(Po - A0V2 -^ai(0)6y 2 X o(^2)(^+3(t2) + (6fc + 6)Vfe+i(i 2 ) 
+ 12fc 2 Vfc-i(i 2 ) + 8fc(fc - l)(fc - 2)^ fc _ 3 (t 2 )) ■ 



We can take 



1 1 

¥>2 =^ai(0)fo o 1/2 Xo(^)(gV-fc+3(t2) + (3fc + 3)Vfc+l(t 2 ) 
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- 6fc 2 ^-i(t 2 ) - h(k - l)(fc - 2)^_ 3 (i 2 )) • 

Step 4 

As fourth equation we obtain: 

(P - Ao)y3 = hfo - P3V0 ■ 
The computation of Pspo gives: 

P^o = - ^a'i(0)&oXo(^)i^fc(*2) 

= - ^K(0)&S /2 XoM(^+3(*2) + (6fc + 6)Vfc+i(t 2 ) 
+ 12fc 2 ^-i(t 2 ) + 8fc(fc - l)(fc - 2)^ fc _ 3 (*2)) • 

Thus, we have 

A3 = 0, 

and 

(P - A )<^ 3 =^m' 1 (0)6y 2 Xo(<T2)(V'fe+3(t2) + (6fc + 6)^+1^2) 
+ 12fcVfc-i(*2) + &k(k - l)(fc - 2)^fc-3(*2)) • 

We can take 

¥>3 =^«a'i(0)foo 1/2 Xo(^2)(^fe+3(i 2 ) + (3fc + 3)^ +1 (t 2 ) 

- 6^-1^2) - ^fc(fc - l)(fc - 2)^ fc _ 3 (t 2 )) • 

Step 5 

The fifth equation reads: 

(Po - ^o)f4 = -Wo - Pa<Po - P2P2 ■ 
The computation of P^ipo gives: 

PWo =1 (fo(0) - 2b (a 2 (0) - §«i(0) 2 )) M^M^XoM 

+ \ (02(0) - ^«i(°) 2 ) V>fc(i 2 )xoM + ^ai(0)t 2 ^fc(*2)ff2Xo(o- 2 ) 

+ ^ai(0)6o ^ 2 A 2 - Vfc(* 2 )XoM + ^i'(0)6^|Vfe(i 2 )Xo(^2) . 

Multiplying by V'fc and integrating with respect to t 2 gives: 

||#~j|2 y P i l Po{ (T 2,t2)lpk{t2)dt2 



2 

Compute 



^/3 2 (0) - 1 (a 2 (0) - | ai (0) 2 ) ) (2/,-' -I- 2k - l) V) (a 2 l 
+ £ f« 2 (0)-|«i(0) 2 ) Xo(<t 2 )- 



P 2 (^ 2 =^ai(0) 2 foo /2 Xo(^)(^ 2 ^+3(i 2 ) + (3fc + 3)tlVfe+i(* 2 ) 
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- 6fc 2 ^Vfc-i(*2) - ^k(k - l)(fc - 2)t^ fe _ 3 (i 2 )) • 
Multiplying by and integrating with respect to t 2 gives: 

/ ^(0*, fc)^) dfe = - (gfc 2 + H fc + 11) a l( 0) 2 X oM • 

We obtain 



A4Xo(ct 2 ) - JiPifo + P2<P2)?Pk dt 2 



= A4X0M - ^/3 2 (0)(2fc 2 + 2fc + l)xoM - (k 2 + k) (a 2 (0) - J«i(0 r ) \ (1 (a 2 ) . 
which holds for any function xo, if we put 

A 4 = ^/3 2 (0)(2fc 2 + 2k + 1) + (fc 2 + fc) (a 2 (0) - lai(0) 2 
To find ip4, we observe that the right hand side has the form 

A 4 </?0 - P4<P0 - P2f2 

= Ao{t2)X(){o2) 77^01 (0)(^fc+i(i 2 ) + 2/c-0 fc _i(t 2 ))cr2Xo(CT 2 ) 

- ^a'i'(0)&o /2 (^ +3 (t2) + (6fc + 6)^ +1 (t 2 ) 
+ 12fc 2 Vfe-i(i 2 ) + Sk(k - l)(fc - 2)Vfe-3(*2))xo(^), 
where A)(*2) has a form 

3 

Ad(*2) = E a m tp k +2m(t2), ao = 0. 

m=— 3 

Thus, we obtain 

¥4 =Bo(t 2 )Xo{o2) - ^72" a l(°) Q^fc+l(*2) - fc^fc-l(<2)) 0-2X0(^2) 



^ a i'(0)^ /2 (^fc+3fe) + (3fc + 3)^ +1 (i 2 ) 



6fc 2 ^-i(t 2 ) - -k(k - l)(k - 2)^_ 3 (i 2 ))xo(^2) , 



3 



4 

where 

B °(^)= E ^^k + 2m{t2). 
m= — 3 

Step 6 

The sixth equation reads: 

(P - A )<^5 = A 5 </? - -P5<A) - ^3^2 - -P2^3 • 

We don't need an explicit formula for (p$ but only its existence. Therefore, we 
only find A5 from an orthogonality condition. We have 

flspo = - h&(0) - 6 (a 2 (0) - lai(O) 2 ) J t 3 2 ip k (t 2 )a 2X o(<T2) 
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23 

^(0)-2& ( 4(0) --a x (0)^(0) 



4^fc(*2)Xo(o-2) 



^ia'('(0)blt 3 Mt 2 )x^'(cr 2 ) + \a'((0)b (2t 3 2 D t2 - 3it 2 2 ) ^fe)xd(^) 



1 



1 



1 



+ f K(0) - -a 1 (0)4(0)J VfcfekoM- 

Multiplying by V>fc and integrating with respect to t2 gives: 
1 



\\Hu\\ 2 



P5<Po(<72,t2)^k(t 2 ) dt 2 

%(0) 



Ah 
11 
32 



^(2fc 2 + 2fe + 1) - ^4(0) - H| ai (o)4(o) ) (fc 2 + fc) 

V 16 / 



oi(0)oi(0)Xo(^2) • 



Next, 



-P3V2 



^ i ai(0)4(0)6^ /2 (^^ fc+3 ( i2 ) + (3fc + 3)t^ fc +i(t 2 ) 
6fc 2 i^fc-ife) - ^fc(fc - - 2)^ fc _ 3 (i 2 ))XoM ■ 



Multiplying by and integrating with respect to t 2 gives: 
1 



ll^ll 
Finally, 



P3 l P2(o- 2 ,t 2 )ipk(t2) dt 2 = - 



15 2 15 11 
32 fc + 32 fc+ 64 



iai(0)4(0)x' (a 2 ) 



-P2V3 



^a 1 (0)a' 1 (0)b 3 /2 q4ip k+ 3(t 2 ) + (3k + 3)t 3 2 ^ k+1 (t 2 ) 
6fc 2 i^fe-i(i 2 ) - ^k(k - l)(k - 2)t 3 2 ip k -3(t 2 ))x'o((T 2 ) . 



Multiplying by V>fc and integrating with respect to t 2 gives: 



1 



\\H k \\ 2 
Using 



P 2 ^3(o- 2 ,t 2 )lpk(t2) dt 2 = - 

, we obtain 



15 2 15 11 

32 + 32 + 64 



mi(0)4(0)x'o(a2) 



p02 J( P W0 +P3V2 +P2<P3)lpk dt 2 

'^-(2k 2 + 2k + 1) - ^4(0) - ioi (0)4(0)^) (fc 2 + fc) 



4fen 



Xo(^) = 0. 



Therefore, the orthogonality condition holds for any function Xoi if we P u t 

A 5 = 0. 

Step 7 

The seventh equation reads: 



(P a - \q)<pg = A 6 <^o - Pefa - P^2 ~ P3V3 - P2<Pi ■ 
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We have 



+ 



#f(0) 1 
3 

3 



66 



« 2 '(o)-f«i(oK(o)-f K(o)) 2 



]a?(0)b ( t* 2 D t2 - ) Vfe(i 2 ) XoM 



+ 



fe(o) 

6o 



i 



- o fl i(°) 2 ) *2^fc(i 2 ) - — a' 1 (0)(2t 2 A 2 - i)^fc(* 2 ) 



o- 2 Xo(o- 2 ) • 



Multiplying by ipk and integrating with respect to t 2 gives: 



+ 



/3 2 '(0) 



(2fc 2 + 2/c + 1) (yfc 2 + fc)a' 2 '(0) 

o 2 

23 



23 2 23 11 
32 + 32^+64 



2bl 



Next, we have 

P 2 ^Pi = 



(2fc + l) 



166 



ai (oK(o ) + ^ + |. + H 

(2fc+l) ai (0) 2 V 2 Xo(^ 2 ). 



K(o)) 2 x£M 



-ai(0)&^ LB (* 2 )xoM 

«i(°) ( ^-0fc+i(t 2 ) - kip k -i(t 2 ) ) (73X0(0-2) 



4^ 2 



^ai'(0) & o /2 (^fc +3 (t 2 ) + (3fc + 3)^ fc+1 (t 2 ) - efc 2 ^!^) 



32 

4 



-fc(fc-l)(fc-2)^ fe _ 3 (t 2 ))xoV 2 ) 



Observe that 



(tifl ,^fc>= 51 £^(4i>k+2m,i>k) = 0. 
m=— 3 

Multiplying by V'fc and integrating with respect to t 2 gives: 



ll^fcll 
Next, 



jr-p y P2¥i{v2,t 2 )^k{t2) dt 2 



32b 



(2fc + l)( ai (0))V 2 X o(a 2 ) 



+ 



«i(oK(o)&o(gfc 2 + gfc + ^)x^'(- 2 ). 



J P3^3=^K(0)) 2 &o /2 i 2 (^V' fe+ 3(t2) + (3fc + 3)^ +1 (f 2 ) 

- 6fc 2 ^-i(i 2 ) - ^(fc - l)(fc - 2)^_ 3 (i 2 ))xo(^ 2 ) • 
Multiplying by ipk and integrating with respect to t 2 gives: 

y P 3 ^(<7 2 ,i 2 )V fe (i 2 ) dt 2 - (|fc 2 + |fc + K(0)) 2 Xo-'M • 
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Finally, the operator P4 has the form 

Pa = Ai(t 2 ) + i«i(0)t 2( r 2 2 - \a'((0)b 2 t s 2 Dl 2 . 

Therefore, 

P W2 =B 4 (i 2 )xo(T2) + ^(ai{0)) 2 b- 1/2 t 2 ^ k+3 (t 2 ) + (3k + 3)Vfc+i(*2) 
- 6fc 2 Vfc-i(*2) - ^(fc - l)(fc - 2)Vfe- 3 (t2))<7|xo(^2) 
+ ^ai(0)ai(0)6g /2 ^(^ fe+3 (t 2 ) + (3fc + 3)^ fe+1 (t 2 ) 



64 
6fc 2 



lk-i(t 2 ) - ^fc(fc - l)(fc - 2)^ fc _ 3 (t2))xo(^) • 



It is easy to see that the term Bi(t 2 ) has the form 

3 

Bi{t 2 )= ^2 7m^fc+2m+l(t 2 ) • 
m=— 4 

Therefore, we have 

(B 4 ,Vfc) =0. 

Multiplying by tpk and integrating with respect to t 2 gives: 
||g fe ||2 y p 4^ 2 (o- 2 ,i 2 )^fc(i 2 )cft 2 

= 32^ (2A: + ^M )) 2 ^^) + ai (o) a '/(o)(gfc 2 + gfc + ^)x£M • 

The orthogonality condition gives 

A 6 Xo(o- 2 ) - ||^||2 J( P 6fo + P4f2 + PWi + Pi<Pi)i>k dt 2 



A 6 Xo(o- 2 ) + ^ 



46 



(2fc 2 + 2fc + 1) 



(A' 2 + / 1 )(a' 2 '(0)-la 1 (0K(0)-i(a' 1 (0)) 2 ) V ;,V 2 ) 



fe(0) 
26 2 



(2A: + l)a 2 xo(a 2 ) 



which has a nontrivial solution xo , if ^6 is an eigenvalue of the operator 



D k = -- 



1 



4r^(2fc 2 +2fc+l)-(fc 2 +fc) (a' 2 '(0) - ioi(0)oi'(0) - J(a'i( )) 



4&o 



+ 



2 

/82(0) 
26 2 



(2fc + l)a| 



Thus, we can take 



26 



®l(2ife 2 + 2*+l) 



- (fc 2 + k) 4'(0) - - fll (0)<(0) - -(ai(0)) : 



1/2 



ft(0) 1/2 (2* + l) 1/2 (2j + l), 

jeN, 
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and 

Xo(fa) = Vjkfa) , 

where "J^fe is the normalized eigenfunction of the operator Dk associated with the 
eigenvalue A 6 . 

Thus, for any j G N, we have constructed an approximate eigenfunction ip^ k of 
the operator H h, ° in the form 

6 

f=o 

such that 

¥?o(*2,cr 2 ) = __^ fc (t 2 )^ ii .( (T2 ) . 
with the corresponding approximate eigenvalue 

A'' fc (fc) = XV /4 Af ■ 

Then we have 

ff h ' ^-A J " fc (/»)^ = O(ft r / 4 ). 

Observe that X ]k (h) = h~ 1 fi jk (h), where fi ]k (h) are given by (|4.3[) . 

The constructed functions <^ fc have sufficient decay properties. Therefore, by 
changing back to the original coordinates and multiplying by a fixed cut-off function, 
we obtain the desired approximate eigenfunctions Uj k , completing the proof of 
Theorem 14.11 

5. Periodic case and spectral gaps 

In this section, we apply the previous results to the problem of existence of gaps 
in the spectrum of a periodic magnetic Schrodinger operator. For related results 
on spectral gaps for periodic magnetic Schrodinger operators, see [8] and references 
therein. 

Let M be a two-dimensional noncompact oriented manifold of dimension n > 2 
equipped with a properly discontinuous action of a finitely generated, discrete group 
r such that M/T is compact. Suppose that H 1 (M, R) = 0, i.e. any closed 1-form 
on M is exact. Let g be a T-invariant Riemannian metric and B a real-valued 
T-invariant closed 2-form on M. Assume that B is exact and choose a real-valued 
1-form A on M such that dA = B. Write B = bdx g , where b £ C°°(M) and dx g is 
the Riemannian volume form. Let 

b = min \b(x)\ . 

Assume that &o > and there exist a (connected) fundamental domain T and a 
constant eo > such that 

\b{x)\ >bo + e , x G dJ- . 

We will consider the magnetic Schrodinger operator H h as an unbounded self- 
adjoint operator in the Hilbert space L 2 (M). Using the results of [7], one can 
immediately derive from Theorem 12.11 the following result on existence of gaps in 
the spectrum of H h in the semiclassical limit. 
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Theorem 5.1. 

Assume the set {x G M : \b(x)\ = bo} contains a smooth curve 7 such that in some 
neighborhood of 7 we have 

C^dix^f < [b(x)[-b Q <Cd(x, 1 ) 2 

with a constant C > 0. Then, for any natural k and N , there exists h^,N > such 
that the spectrum of H h in the interval 

[(2k + l)hb Q + h 2 m k , (2k + l)hb + h 2 M k ] , 

where [m k ,M k ] is the range of the function Vk on 7 defined by (|4.1j) . has at least 
N gaps for any h g]0, hk,N\ ■ 

Similarly, Theorem 14.11 implies the following result. 
Theorem 5.2. 

Assume the set {x £ M : \b(x)\ — bo} contains a smooth curve 7 such that in some 
neighborhood of 7 we have 

C- 1 ^^) 2 < \b(x)\-b <Cd(x,j) 2 

with a constant C > 0, and for some k € N, there exists a unique minimum xq G 7 
of the function Vk on 7 defined by (|4.1I) . which is nondegenerate : 

5 k = V k "(x )>0. 

Then, for any natural N, there exist 

Ck,N > ^-6l /2 ^ a /2 (2k + l) 1 / 2 (2N + 3) 
and h k ,N > such that the spectrum of H h in the interval 

[(2k + l)hb + h 2 m k + h^ 2 ^-S 1 k /2 f i 1 /2 (2k + I) 1 ' 2 , (2k + l)hb + h 2 m k + h 5 / 2 C k . N ] 
has at least N gaps for any h G]0, h k ^] ■ 

Appendix A. Some facts from geometry 

Here we refer to [B] for more material. Let M be a Riemannian manifold. Denote 
by V the Levi-Civita connection associated with the Riemannian metric g. Recall 
that a connection in the tangent bundle of M is a map 

V : C°°(M,TM) x C°°(M,TM) -» C°°(M,TM), (X,Y) h-> V x Y , 

which satisfies 

(1) for any X € C°°(M,TM) the map V x : C°°(M,TM) C°°(M,TM) is 
linear; 

(2) for any Xi,X 2 £ C°°(M, TM) and Y £ C°°(M, TM) 

V Xl+ x 2 s = V Xl Y + Vx 2 Y:, 

(3) for any / G C°°(M), X G C°°(M,TM) and Y G C°°(M,TM) 

V fx Y = fV x Y; 

(4) for any / G C°°(M), X G C°°(M,TM) and Y G C°°(M,TM) 

V x (fY) = fV x Y + X(f)Y . 
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The Levi-Civita connection is the unique connection V in TM, which is compat- 
ible with the Ricmannian metric (with the corresponding scalar product denoted 
by (-,-)) in the sense that, for any X,Y,Ze C°°(M, TM) 

X(Y,Z) = (VxY,Z) + (Y,V x Z), 

and torsion-free, in the sense that, for any X, Y e C°°(M, TM) , 

VxY -V Y X - [X,Y] = 0. 

The curvature of the connection V is the operator 

R(X,Y) : C°°(M,TM) ->■ C°°(M,TM) , 

associated with any pair X, Y e C°°(M,TM) by the formula 

R{X, Y)Z = V X V Y Z - V Y V X Z - V [X , Y ]Z 

for any Z G C°°(M,TM). 

In a local coordinate system (x 1 , x 2 , . . . , x"), one can write 

where 

ftri fir* ™ ™ 
£j7£ ~ ~dxT ~dx k ~ ^ u jm Z^i it km ■ 

m—l m—l 

The Ricci tensor of M is given by 

n 

RiCafc = Rltb ■ 

i=l 

The scalar curvature of M is defined by 

K = ^.g afc Ric ah . 

ab 

The Riemannian curvature tensor is a R of type (4,0), defined in local coordinates 
as 

n 

Rijki = ^ 9imRJk£ ■ 
m—l 

One can also give its invariant definition. The corresponding map 
R : X(M) x X(M) x X(M) x X(M) -> C°°(M) 

is given by 

R(X U X 2 ,X 3 ,X 4 ) = {RiX^X^X^Xi), X 1 ,X 2 ,X 3 ,X 4 e X(M). 

If dim M = 2, then the Riemannian curvature tensor R has 4 nontrivial compo- 
nents 

-R1212 = —R2W2 = R\221 = — -Rl221i 
Other components equal zero. We have 

2i?1212 = R(gn922 — ffi 2 )) 
moreover the scalar curvature k is related with the Gauss curvature K by 

R = 2K. 
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Now assume that M is two-dimensional, and 7 is a one dimensional smooth 
submanifold (a closed curve). Let s be the natural parameter along 7. Assume 
that 7 is oriented, and choose the external unit normal vector TV at each point of 7. 
In a tubular neighborhood U of 7, consider the associated Fermi coordinate system 
(Xq,Xi) = (s,t), which is defined as follows. For any x € U, t is the distance from 
x to 7 and s is the coordinate of the intersection point of the minimal geodesic £, 
passing through x orthogonally to 7. Then at each point x e U the vector A = 
coincides with the unit tangent vector of the minimal geodesic £, passing through 
x orthogonally to 7. Thus, by definition, we have 

||A|| = 1, V N N = Q. 

Consider the vector field AmU defined in local coordinates (X ,Xi) = (s,t) as 




For any x £ 7 the vector A(x) is tangent to 7 and 

\\A(x)\\=l. 

Observe that 

[A,N}=0. 

Then we have 

V N V N A + R AN N = 0. 
Indeed, by definition, we have 

R AN N = Va^nN - Vn^aN - V [A , N] = -V N V A N = -V N V N A. 
Next, we have 

N(A, N) = (W N A, N) + (A, Vat TV) 
= (V N A,N). 

On the other hand, since ||A|| = 1, we obtain 

0=A(N,N) =2(V A N,N). 

Therefore, we have 

N(A,N)=0, 

that is, (A, N) is constant along the integral curves of A, which are minimal 
geodesies orthogonal to 7. Since (A, A) = on 7, it follows that 

(A,N)=0. 

Thus, the coefficients of the metric have the form 

g 00 (s,t) = (A,A):=a(s,t), g 01 (s,t) = (A,N) = 0, g n (s,t) = (A, A) = 1. 
Observing that 

a(s,0) = 1, 

the metric g has the form 

g — a(s, t)ds 2 + dt 2 . 
Let us compute Taylor's expansion of a at t = 0: 

g Q0 (s, t) = l + ai(*)t + a 2 {s)t 2 + 0(t 3 ). 
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First, we recall the definition of the mean curvature of 7. Since \\A\\ — 1 on 7, we 
have 

k = (T a A,N) = {V A A,N). 

Thus we obtain that 

(A.l) ai(s) = N(A,A) = 2{V N A,A) = 2(V A N,A) = -2(V A A,N) = -2k. 
Then, 

02(5) =\n 2 {A, A) = (VatVatA, A) + \\V N A\\ 2 . 

For the first term, we get 

{V N V N A,A) = -(R AN N,A) = -R 1001 = ~R 

For the second term, since (^4, N) = on 7, we obtain 

||VatA|| 2 = (W N A,A) 2 + (W N A,N) 2 = k 2 . 

Thus, we have 

(A.2) a 2 = -^R + K 2 . 
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